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ABSTRACT 



Hydrogen atoms inside virialized minihalos (with Tvir ^ lO^K) generate a radia- 
(N ■ tion background from redshifted 21-cm line emission whose angular fluctuations reflect 

clustering during before and during reionization. We have shown elsewhere that this 
emission may be detectable with the planned Low Frequency Array (LOFAR) and 
^SJ ' Square Kilometer Array (SKA) in a flat Cold Dark Matter Universe with a cosmo- 

0^ \ logical constant (ACDM). This is a direct probe of structure during the "Dark Ages" 

, at redshifts z ^ 6 and down to smaller scales than have previously been constrained. 

' In our original calculation, we used a standard approximation known as the "linear 

bias" [e.g. Mo & White (1996)]. Here we improve upon that treatment by considering 
the effect of nonlinear clustering. To accomplish this, we develop a new analytical 
method for calculating the nonlinear Eulerian bias of halos, which should be useful for 



O ' other applications as well. Predictions of this method are compared with the results of 

^ , ACDM N-body simulations, showing significantly better agreement than the standard 

^ ' linear bias approximation. When applied to the 21-cm background from minihalos, our 

formalism predicts fluctuations that differ from our original predictions by up to 30% 
at low frequencies (high-z) and small scales. However, within the range of frequencies 
^ , and angular scales at which the signal could be observable by LOFAR and SKA as 

' currently planned, the differences are small and our original predictions prove robust. 



Our results indicate that while a smaller frequency bandwidth of observation leads to 
a higher signal that is more sensitive to nonlinear effects, this effect is counteracted by 
the lowered sensitivity of the radio arrays. We calculate the best frequency bandwidth 
for these observations to be Ai'obs ~ 2 MHz. Finally we combine our simulations with 
our previous calculations of the 21-cm emission from individual minihalos to construct 
illustrative radio maps at z = 9. 

Key words: cosmology: theory — diffuse radiation — intergalactic medium — large- 
scale structure of universe — galaxies: formation — radio lines: galaxies 



1 INTRODUCTION 

Despite striking progress over the past decade, cosmologists have still failed to see into the "Dark Ages" of cosmic time. 
Although the details of this epoch between recombination at redshift z ~ 10^ and reionization at 2; > 6 are crucial for 
understanding issues ranging from early structure formation to the process of reionization itself, no direct observations of 
any kind have been made in this redshift range. Recently, we made the first proposal for such direct observations (Iliev et al. 
2002, from hereafter Paper I) based on coUisional excitation of the hydrogen 21-cm line in the warm, dense, neutral gas in 
virialized minihalos (halos with virial temperature Tvir ^ 10*K), the first baryonic structures to emerge in the standard CDM 
universe. We showed that collisional excitation is sufficient to increase the spin temperature of hydrogen atoms inside minihalos 
above the temperature the cosmic microwave background (CMB). Minihalos should generally appear in emission with respect 
to the CMB, p roducing a background "21-cm forest" of redshifted emission lines, well-separated in frequency. Unlike all 
previous works jHoean fc Reedl 19791: IScott fc Ree3ll990l : ISubramanian fc Padmanabhanlll993t iMadau. Meiksen fc Ree3ll997l : 
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IShaver et alJll999l : iTozzi et al.ll2000D , this mechanism does not require sources of Lya radiation for "pumping" the 21-cm hne 
and decouphng it from the CMB. 

In Paper I we calculated the 21-cm emission properties of individual minihalos in detail, along with the total background 
and its large-scale fluctuations due to clustering. Although individual lines and the overall background are too weak to be 
readily detected, the background fluctuations should be measurable on ~ 10' — 100' scales with the currently-planned radio 
arrays LOFAR and SKA. We demonstrated that such observations can be used to probe the details of reionization as well as 
measure the power spectrum of density fluctuations at far smaller scales than have been constrained previously. 

In the current paper we extend these results to smaller scales. Our previous calculations showed that the fluctuation 
signal increases as the beam size and frequency bandwidth decrease, which corresponds to sampling the 21-cm emission from 
minihalos within smaller volumes. As these volumes correspond in turn to length scales that are more nonlinear, nonlinear 
effects have the greatest impact on the angular scales and frequency bandwidths at which the signal is the strongest. 

This issu e is of particular importance as our previous investigations relied on the standard, simplified nonlinear bias of 
IMo fc Whit j il996l) . which breaks down at many of the relevant redshifts and length scales. For example, an rms density 
fluctuation in a Cold Dark Matter universe with a cosmological constant (ACDM) at z — 6(8), is strongly nonlinear {a{AI) ^ 1) 
for M ^ W^Mq (6.2 x IO'^Mq), which roughly corresponds to the region sampled by beam sizes 9" (3"), for frequency 
bandwidths 12 kHz (3 kHz), respectively. Nonlinear effects can influence the predicted background fluctuations on even larger 
scales, up to few comoving Mpc (corresponding to few hundred kHz frequency bandwidths and 1-10 arc min beams), even if 
these scales are still not strongly nonlinear. Additionally, rare halos are always more strongly clustered than the underlying 
density distribution (i.e. the bias is > 1), again bringing nonlinear issues to the fore. 

In order to address these important issues we have carried out a series of high-redshift N-body simulations of small scale 
structure formation, which we use to construct 21-cm line radio maps that illustrate the expected fluctuations in the emission. 
As no current simulations are able to span the full dynamic range relevant to 21-cm emission, however, we extend our results 
by developing a new formalism for calculating the nonlinear Eulerian bias of halos, which is based on the Lagrangian bias 
formalism of Scannapieco & Barkana (2002, hereafter SBQ2). We describe this approach in detail in this paper, verify it by 
comparing it with the results of our N-body simulations, and apply it to to calculate improved predictions for the fluctuations 
in the 21-cm emission. These are then compared to the predictions given in Paper I, quantifying the impact of nonlinear 
effects on minihalo emission. 

The structure of this work is as follows. In §|5|we describe a set of simulations of minihalo emission at high redshift, and 
use these to construct simulated maps at small angular scales. In §|31we present our improved calculation of the Eulerian bias 
and verify it by comparison with the numerical simulations presented in §|21 In §0]we modify the calculation of the radiation 
background from minihalos to incorporate the contribution due to nonlinear clustering of sources, according to the formalism 
described in § |3 and describe the results of our nonlinear formalism. Conclusions are given in § |^ 



2 NUMERICAL SIMULATIONS AND SIMULATED 21-CM RADIO MAPS 

We simulated the formation of minihalos in three cubic computational volumes of comoving size 1 Mpc, 0.5 Mpc, and 0.25 Mpc, 
respectively. We used a standard Particle- Particle/Particle- Mesh (P'^M) algorithm (Hockney & Eastwood 1981), with 128"^ 
equal-mass particles, a 256^ PM grid, and a softening length of 0.3 grid spacing. Here and throughout this paper we consider a 
flat ACDM model with density parameter Qo — 0.3, cosmological constant Aq — 0.7, Hubble constant Ho = 70 kms~^Mpc~^, 
baryon density parameter Q,i, — 0.02h~^ (where h = _H'o/100kms~'^Mpc~^), and no tilt. The initial conditions where generated 
using the transfer function of Bardeen et al. (1986) with the normalization of Bunn & White (1997) (for details, see Martel 
& Matzner 2000, §2). All simulations started at redshift 2 = 50 and terminated at redshift z = 9. To identify minihalos, we 
used a standard friends-of-friends algorithm with a linking length equal to 0.25 times the mean particle spacing. We rejected 
halos composed of 20 particles or less. Table 1 lists the comoving size Lbox of the box, the comoving value of the softening 
length rj, the total mass Aftot inside the box, and the particle mass A/part. 

Once the halos are identified, we can use their properties and distribution to compute the corresponding radio maps. We 
assign to each halo a 21-cm fiux based on its mass and redshift of formation by modeling the halos as Truncated Isothermal 
Spheres (Shapiro. Ilicv & Raea 1999; Iliev & Shaniro 2001}. The 21-cm fluxes from individual minihalos are obtained by 
solving the radiative transfer equation self-consistently through each halo to obtain the line-integrated flux, as described in 
Paper I. Let us consider a cylindrical volume with radius R = A6'beam(l + z)Da{z)/2, and length Lbox , where Da{z) is 
the angular diameter distance, A^boam is the beam size of the observation. Then the beam-averaged differential antenna 
temperature STt is calculated using equation 6 of Paper I, which can be discretized as follows; 



where H{z) is the Hubble constant at redshift z, uq is the rest-frame line frequency, 7r_R^Lbox is the comoving volume of 




(1) 
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Table 1. Parameters of the N-body simulations 



ibox (Mpc) 


V (kpc) 


Mtot (Mq) 


Afpart(M0) 


1.0 


1.172 


4.079 X lO^" 


1.945 X 10" 


0.5 


0.586 


5.099 X 10^ 


2.431 X 10^ 


0.25 


0.293 


6.374 X 10** 


3.039 X 10^ 



the beam (Ai/ofr)i, {STb,,^g)i, and Ai are the effective Une width, the line-center differential brightness temperature, and the 
geometric cross-section of halo i, respectively, and we use Gaussian filter in order to ensure that the contribution of a particular 
minihalo to a particular pixel varies smoothly with the location of the minihalo. 

The resulting maps are shown in Figure Due to the small box sizes that are required to resolve the minihalos, the 
beams used to produce the maps are also very small, ranging from 4" to 0.25". The fiuxes from such small beam sizes are well 
below the sensitivity limits of the currently planned radio arrays LOFAR and SKA. Additionally, the larger-box (1 Mpc and 
0.5 Mpc) simulations do not have sufficient mass resolution to resolve the smallest halos, while the larger-mass minihalos are 
not present in the smaller box simulations. Therefore the fiux levels are somewhat underestimated in all cases, and should be 
considered only as illustrative of fluctuations on very small scales. 



3 NONLINEAR CLUSTERING REGIME AND COMPARISON TO SIMULATIONS 
3.1 An analytic approach to nonlinear clustering 

As our simulations are unable to span a sufficient range of scales to both resolve all minihalos and reproduce the physical 
scales accessible to observations, we adopt instead an alternate technique to estimate the full nonlinear signal. We apply 
an approximate formalism that is an extension of the one developed in SB02. A more complete discussion of this method, 
refinement of the basic formalism, and detailed comparisons with simulations will be presented in a separate paper. In this 
section we summarize the basic features of this approach and show th at it is sufficient for our purposes in this investigation. 
For an alternative approach see lCatelan. Matarese. fc Porcianil ( 199jj^ . 

In SB02 the authors extended the Press-Schechter formalism as reinterpreted by iBond et alJ (Il99lh to calculate the 
number density of collapsed objects in two regions of space initially separated by a fixed comoving distance. In this approach 
virialized halos are associated with linear density peaks that fall above a critical value usually taken to be 5c = 1.686. Here 
the linear overdensity field is defined as 5(x)D~^(z) = p{:x., z) / p{z) — 1, where p(x) is the linear mass density as a function 
of position and and p is the mean mass density. Note that here the linear overdensity field 5(x) is described in terms of its 
value extrapolated to z = 0, and its evolution with time is subsumed by the "linear growth fact or" D(z) = 5(0)/S(z). Other 
choices for 6c and i ts evolution have been explored in the literature, by fitting to simulations (|S heth^ Mo^ fc Tormenl EoOll : 
Ijenkins et al. lEoOlh or incorporating the Zel'dovich or adhesion approximations ihee fc Shandarinl^99nMencill20MF !^While 
these methods improve the accuracy of the Press-Schechter technique, they provide no information about nonlinear clustering. 
Here we strive to develop a formalism that is applicable to nonlinear clustering regardless of the recipe chosen for Sc. 

Once 5c determined, the presence of a halo can be associated with a random walk procedure. At a given redshift we 
consider the smoothed density in a region around a point in space. We begin by averaging over a large mass scale M, or 
equivalently, by including only small comoving wavenumbers k. We then lower M, adding k modes, and adjusting 5{M) 
accordingly. This amounts to a random walk in which the "time variable" is the variance associated with the filter mass and 
the "spatial variable" is the overdensity itself. This walk continues until the averaged overdensity crosses 5cD-^{z) and we 
assume that the point belongs to a halo with a mass M corresponding to this filter scale. 

The problem of halo collapse at two positions in space can be similarly associated with two random walks in the presence 
of a barrier of height 5cD~^{z). These random walks are extremely correlated at large scales, and become less and less so as 
the modes of smaller and smaller scales are added. In SB02 this process was approximated by a completely correlated random 
walk down to an overall variance 5* = ^, where ^ is the correlation function that accounts for the excess probability of finding 
a second halo at a fixed distance from a given halo. This correlated random walk was then followed by uncorrelated walks 
down to the final variances associated with the masses of the objects, S = a^{Mi) and a^{M2). 

This procedure lends itself directly to calculating the overall increase in number density due to infall. The fact that the 
overdensity on large scales will contract the distances between halos provides an additional contribution at the length scale 
X such that o-^{x) = ^, the scale down to which the evolution of the two points is completely correlated. Let 5o be the linear 
overdensity at redshift zero at that scale. The probability distribution of 5o, Qo is then simply given by a Gaussian, with a 
barrier imposed at ^{z) = 5cD~^{z): 

Qo{i^,5o,0 = [G{5o,0 ~ G{2u - 5o,O]0{u - So), (2) 
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Figure 1. (upper panels) Differential brightness temperature radio map of the 21-cm emission from minihalos at 2 = 9, produced as 
described in the text from ACDM N-body simulation of: (upper left panel) box size 1 Mpc and beam size A^bcam = 4", and (upper right 
panel) box size 0.5 Mpc and beam size A^bcam = 2". (lower left panel) box size 0.25 Mpc and A^bcam = l"i '^'nd (lower right panel) 
box size 0.25 Mpc and Aeheam = 0.25". 



Nonlinear clustering of minihalos 5 



where G{S,^) = {2n^)^^^^e^^ , and 9 is the Heaviside step function. 

We now consider f2{u{z), fi, af, 5o, daf dal, the joint probability of having point A in a halo with mass corresponding 
to the range af to + dai and point B in a halo with mass corresponding to the range erf to o"! + da2, whose random walks 
pass though a value of So at the a'^{M) — ^ scale. In this case we obtain an expression analogous to equation (37) of SB02 



dSi 



dSi 



d52 



d52 



)i2(i^, 5i,52,o-l,al,i{r)), 



(3) 



where now Qr2 = G{5i — 5o,af — ^)G{52 — So, o"! ~ C)i ^-nd the reader is referred to SB02 for a more detailed derivation of /2 
from the underlying probability distribution. 

Averaging equation JSJ over the probability distribution for So as given in equation ((SJ and carrying out the partial 
derivatives we obtain the Eulerian "bivariate" number density, the joint probability of having point A lie in a halo in the mass 
range Mi to Mi + dMi and point B at a comoving distance r lie in a halo in the mass range M2 to M2 + dA/2 at a redshift z: 



'2„2 



12, E 



(r,z) 



dMidM2 
where 



P 

Ml 



dal 



dMi 



M2 



dcTo 



dM2 



h,E{v{z),Gi,a2,i(r)), 



dSoQo{5o,0 g{D~'So) f{v - 5o, e - ^l) g{D-'So)f{iy - <5o, C - '^1) 



(4) 



(5) 



and f{u,a') = {2-.y^/\u / a') 



' /2a 



Here g{D ^So) is p/p at the scale at which the points A and B are completely 



correlated. This can be thought of as the contraction of a large spherical perturbation that surrounds both points and 
contains a total mass M of material, where (j^(M) = ^. This contribution is then just (1 + D~^So) in the linear regime, and 
for our purposes here we assume g{D~^So) = (1 + D~^So) for all values as this reproduces well the results from simulations. 
Finally we define the Eulerian bias as 

d^n\2,E f dm^E dn2,E \~^ ^ 

dMidM2 V dMi dM2 ) ^ ' 

where we divide by duE/dM, a single point number density that accounts for the self correlations between the collapsed peaks 
and the overdense sphere, which are over-counted in equation @. To compute this probability we again carry out an average 
over the probability distribution Q, but in this case we consider fi{i>, So, ai)dai — f{v — So,<J^ — £,)dai, the probability of 
having a single point in a halo with a mass corresponding to the range to af + da\ whose random walk passes through So 
at the o-^(M) = ^ scale. This gives duE/dM = (p/M) \da'^/dM\ (1 + d^/dij'^)fE{v,o'^) where 



fE{v,o'^,£) 



dSoQoiSo,i^,OgiSoD ^)f{v - So,i- cr^). 



Combining these expressions yields our final expression for the bias 
f2,E{v, o-f,ai,^) 



1. 



fE{v,(Tl,C) fE[v2,(Tl,C) 

Note that at large distances we can work to order So to determine the asymptotic limit as ^/cr^ 

l-D-Uo ' 



0. In this limit 



g{SoD-^)f{u-So,c7^ -0 



(7) 



(8) 



(9) 



and the only surviving term in equation Q is the cross term between terms of order So- All other terms cancel out between 
the numerator and the denominator, giving 



1 + 



1 + 



Thus our formalism reproduces the usual bias as in IMo fc Whitl il996l) at large distance, as was used in Paper I. 
Finally, we define the flux-averaged correlation between objects as 

/ dMi J dM2 F{Mi){dn/dMi) F{M2){dn/dM2) b%^D-'^ 



[/ dM F{M)(dn/dM)\ 
where F{M) is the line-integrated flux of a minihalo of mass M. 



(10) 



(11) 



3.2 Nonlinear properties and comparisons with simulations 

In Fig. |5|we plot b%^D~^ as calculated from equation and the st andard expression, equation IIUII . using the fit to the 
density fluctuation power spectrum given bv lEisenstein fc HrJ (1 19991 '). In each panel we show the behavior of IO'^Mq and 
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Figure 2. Comparison of standard and nonlinear bias. In all panels, the solid curves give fe|,D~^5 as calculated in equation JSJ while 
the dotted lines give the standard bias expressions as used in Paper I. For each case, the upper curves assume a fixed mass of 10 '^Mq, 
corresponding to a spherical perturbation with a comoving radius of approximately 38 kpc, while the lower curves assume a fixed mass 
of 10^ Mq, corresponding to a comoving radius of 18 kpc. 



IO^Mq objects at four different redshifts, which span the detectable range. At large distances, the clustering expressions 
approach each other asymptotically, as expected from equation (IIUH . On the other hand, the solid lines shoot up dramatically 
at the smallest distances. This is because the separation between the halos is comparable to their radii, and thus the likelihood 
of finding a second collapsed halo at the same point becomes infinite as r — > 0. At intermediate distances the behavior is more 
complex. For relatively rare objects, corresponding to high redshifts, the nonlinear values consistently exceed the standard 
ones. For example, at a distance of 0.1 comoving Mpc, the nonlinear lO^M© value is almost twice that of the standard result 
a.t z — 20, while it is only about 70% of the standard result at z; = 6. 

To compare this behavior to numerical results, we calculated the correlations between halos for each of our simulations. 
In this case b%^D~^ is directly comparable to the halo correlation function ^i,j{r), the excess probability of finding a particle 
within a mass bin rrii and a particle within a mass bin rrij separated by a distance r. This expression is symmetric, so that 
~ Cj, »('■)• We define {Nij(r)) as the mean number of particles in a mass bin rrij, located within a radius r of a particle 
in a mass bin rrii, averaged over all particles in a mass bin A4i. This is is obtained by dividing the total number of pairs with 
separation r' < r, by the number Ni of i particles. 



+ 3 / ^.Ar')r"dr' 



(12) 



where Nj is the number of particles of mass rrij and Vbox is the volume. Finally, we differentiate equation (1121 to get ^i,j{r). 

Using these definitions we examined the distribution of halos in each of our simulations at 2 = 15 and 2 = 9. In each case 
we divided the halos into mass bins containing roughly the same number (~ 600 ) of objects and compared their correlations 
to those given by equation ||HJ. In the left panels of Fig.j^we compare the mass averaged correlation between halos, 

C.. = I / TTiz / ^bl^D-' I I / / -I- I , (13) 





to the numerical value, ^ij, for objects in three mass bins in our 1 Mpc'' simulation at 2 = 15 (Fig. |3J. For all mass ranges, 
there is a good match between the numerical and our analytic expression, with both values equaling or exceeding the 
linear predictions at most distances. The only case in which the predictions of the simulations fall below equation 1101 . is 
at large distances, where the separations between halos begin to approach the simulation box size, and some damping is to 
be expected. Comparisons with our smaller box simulations yielded comparable results, but with larger damping at these 
distances, indicating that the analytic and simulated correlations are a good match at all reliable separations. On the other 
hand, the cross-correlation functions are a poorer match, particularly if the mass bins are very different. For all correlations 
however, our nonlinear estimate falls between the standard and simulated results, indicating that the correlations between 
halos exceed those predicted by equation l|l()|l . and that our approach gives a conservative estimate of this excess. 

In the right panels of Fig.|3 we repeat this comparison with our 1 Mpc'^ simulation at 2 = 9. In this case, we divided our 
halos into five mass bins, three of which are shown in this figure. In this case there is good agreement between the correlations 
and cross- correlations obtained from the simulations and our nonlinear estimates, for all mass scales considered. In fact the 
analytical estimates are consistent with the results of the simulations for all distances that are small relative to the box, apart 
from the dip at small distances in the mi mi case, which is most likely due to small- number statistics. Note that unlike the 
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1.0 Mpc z=15 10 Mpc z = 9 




x/Mpc x/Mpc 

Figure 3. Comparison between nonlinear analytic expressions for the bias and simulations at two different redshifts. In the left panels 
2 = 15 and the halos are divided into three mass bins: 4.1 X IO^Mq ^ mi ^ 5.5 X lO^M©, 5.6 X IO^Mq ^ m2 ^ 1.0 X IO^Mq, and 
1.1 X 10® -Mq ^ m3 7.7 X W'^Mq. In the right panels z = 9 and halos are divided into five mass bins, three of which are shown: 
4.1 X lO^M© ^ mi ^ 4.7 X IO'^Mq, 6.2 X IO^Mq ^ ma ^ 8.2 x W^Mq, and 1.3 X 10*^ Mq ^ ms ^ 3.4 X W^Mq. In both cases, the points 
are the correlation function as calculated from the simulations, while the solid lines give our analytic estimates as given by equation ISt . 
The error bars are 2-a estimates of the statistical noise given by 2^i^i{r)/ ^ Ni^i{< r). Finally, the dashed lines give correlations as 
estimated from the standard approach, equation 1101 . 

higher redshift case, our nonlinear estimates now fall below the standard analytic approach at larger distances, as we saw in 
the lower redshift cases in Fig. |5| 



4 FLUCTUATIONS OF THE 21-CM EMISSION FROM MINIHALOS: EFFECTS OF THE 
NONLINEAR BIAS 

The amplitude of q-a (i.e. q times the rms value) angular fluctuations in the differential brightness temperature STt (or, 
equivalently, of the flux) in the linear regime are given by 

{ST^y^^ = qb{z)apSJl (14) 
(Paper I), where 6(2) is the mean flux- weighted bias, and 

where P{k) is the linear power spectrum at 2 = 0, and the factor (1 + fx^Y, wher e f ~ [r2(2)] °''^ is the correction to 
the cylinder length for the departure from Hubble expansion due to peculiar velocities l)Kaiseilll987h . In order to apply the 
formalism developed in we use that the power spectrum is the Fourier transform of the correlation function ^(r), obtaining 



2 _ 32D-^(2 ) / ^ 2 



drr^ar)f{r,R,L), (16) 



where 
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Figure 4. Predicted 3-cr differential antenna temperature fluctuations at 2 = 7 and z = 8.5 vs. observer frequency bandwidth Ai^obs for 
standard ACDM model using Mo & White bias (long-dashed curves) at angular scales (a) (left) A9 = 9', and (b) (right) A9 = 25'. Also 
indicated is the predicted sensitivity of LOFAR and SKA radio arrays for integration times 100 h (dashed) and 1000 h (dotted) (for right 
panels sensitivity curves for LOFAR and SKA are identical) assuming rms sensitivity oc (see|http:/ /www.lofar.org/science |. 




Figure 5. Predicted 3-cr differential antenna temperature fluctuations at z = 8.5 and z = 20 vs. angular scale Afibeam for standard 
ACDM model for observer frequency bandwidth Au^hs = 2 MHz : using Mo & White bias (long-dashed curves) and current results 
(solid curves). Also indicated is the predicted sensitivity of LOFAR and SKA integration times of 100 h (short-dashed lines) and 1000 
h (dotted lines), with compact sub aperture (horizontal lines) and extended configuration needed to achieve higher resolution (diagonal 
lines) (see http://www.lofar.org/science and Paper I for details). 



Using equations l|14^ and Hlfifl the mean squared angular fluctuations become 

(Sn'} = (^-^y J\rr' [D-''b\z)ar)m?] f{r,R,L). (18) 

We apply the formalism developed in JSlby simply replacing \JD~'^b'^{z)^^{r){5Tbf'^^ with the corresponding quantity D~^b'^^£_{r){SThY 
for the mean Eulerian flux-weighted bias calculated in JHl As we showed, at large distances the two expressions are equivalent, 
while at smaller distances the new formalism better reproduces the numerical results. 

We start with deriving the optimal frequency bandwidth for observing the 21-cm emission from minihalos. As long as 
Ai^obs and the beam size A^bcam are large enough to provide a fair sample of the halo distribution, the mean 21-cm flux is 



Nonlinear clustering of minihalos 9 




Figure 6. (a) (left panels) Predicted 3-cr differential antenna temperature fluctuations at 2 = 8.5 and 2 = 20 at small angular scales vs. 
observer frequency bandwidth Au^ha for standard ACDM model and for angular scale A6 = 1', using Mo & White bias (long-dashed 
curves) and current results (solid curves), (b) (right panels) Predicted 3-cr differential antenna temperature fluctuations at A^bcam = 25' 
vs. redshift 2 for ACDM model, same notation as in a). The predicted sensitivity for integration times 100 h (dashed) and 1000 h (dotted) 
of both LOFAR ("L") and SKA ("S"), are plotted (for bottom panel, sensitivity curves for LOFAR and SKA are identical) assuming 
rms sensitivity oc (see http://www.lofar.org/sciencel. 



independent of Ai^obs, while the fluctuations grow as Ai/obs decreases. However the sensitivity of the radio array deteriorates 
for smaller bandwidths as Ai^^^^ (e.g. Shaver et al. 1999; Tozzi et al. 2000). The results for 3 — u fluctuations at redshifts 
z = 7 and z = 8.5 vs. observed frequency bandwidth A/^obs for beam sizes A^beam = 9' and 25' are shown in Fig. 21 along 
with the results obtained using the linear bias in order to facilitate comparison. As the bandwidth increases, the integration 
time required to detect the signal decreases, reaching minimum at about A/^obs = 2 MHz and remaining roughly unchanged 
thereafter. The fluctuations decrease with increasing bandwidth, however, hence the optimal bandwidth is about Ai/obs ~ 2 
MHz. We see that at large scales the current results are largely indistinguishable from the linear results, as expected, since 
at such large scales the bias calculated in §|21 reduces to the linear expression in equation 1101 . For small values of A^obs the 
differences between the two results grow to ~ 10%. However, in the observable range the two results hardly differ. 

Illustrative results for 3 — a fluctuations at redshifts z = 8.5 and z — 20 vs. the beam size A^bcam for fixed observed 
frequency bandwidth A/^obs = 2 MHz are shown in Fig.|K| The approximate sensitivities for the LOFAR and SKA radio arrays 
are also shown where appropriate (see Paper I for details), but from here on we concentrate on the differences between the 
results from the two approaches and the robustness of our original predictions rather than observability, which was discussed 
in Paper I. At redshifts 2 < 15 using the linear bias gives an overestimate, although a small one compared to the various 
uncertainties of the calculation, while at higher redshifts using the linear bias gives a small underestimate of the fluctuations. 

In Fig. |S| (left panels) we plot the 3 — a fluctuations at redshifts z — 8.5 and z — 20 vs. the observer frequency bandwidth 
Ai/obs for fixed beam size A^bcam ~ 1'. We have chosen such a small beam size in order to investigate the upper limit of the 
deviation of the results due to the nonlinear clustering of halos. For small values of both Ai/obs and AOteam the differences 
are as large as 20-30%. In Fig. |n] (right panels) we plot the 3 — a fluctuations for beam sizes ASbcam = 9' and 25' vs. redshift 
z (i. e. the spectrum of fluctuations) for observer frequency bandwidth Ai/obs ~ 2 MHz. Again, the same patterns emerge, 
where the linear bias approach overestimates the fluctuations at the lower end of the range of considered redshifts by up to 
10 %, while underestimating the fluctuations at the higher redshifts, although by a smaller fraction. 



5 CONCLUSIONS 

We have developed a new analytical method for estimating the nonlinear bias of halos and used it to calculate the fluctuations 
of the 21-cm emission from the clustering of high-z minihalos. This method is likely to be useful in tackling a much wider 
range of problems that lie beyond the capabilities of current numerical simulations, and will be reflned and further verified 
in a future publicat ion. The minihalo bias predicted by our method at large scales reproduces the linear bias, as derived in 
IMo fc Whit j il996l) . and is much larger at small scales, confirming naive expectations. At intermediate scales, however, its 
behavior is more complex and both mass- and redshift-dependent. For very rare halos at high- 2: (roughly z > 15), the standard 
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linear bias is lower than our nonlinear prediction at all length scales. At the lower end of the redshift range we consider (down 
to 2; = 6) the linear bias is higher at intermediate scales (few hundred kpc to few Mpc comoving) and lower at smaller scales. 

We have also compared the predictions of our new method with the results of N-body numerical simulations, which we 
used both to verify our approach and to produce sample radio maps. Due to the limited dynamical range of our simulations, 
however, these maps are only illustrative of the behavior of the fluctuations on very small scales, below the sensitivity limits 
of LOFAR and SKA. On the scales at which the simulations are reliable, we find excellent agreement between these results 
and our analytical approach. The most significant discrepancies occur in the calculation of cross-correlation functions of very 
different mass bins at higher redshifts. However, these departures are relatively modest, and in all cases our method reproduces 
the simulation results significantly better than the linear bias. 

Despite these differences, we find our original linear bias predictions for the fluctuations of the 21-cm emission from 
minihalos to be robust at the scales and frequencies corresponding to observable signals. The prediction using the flux- 
weighted nonlinear bias never departs from the linear prediction by more than few percent in that range, well within the other 
uncertainties of the calculation. This robustness is partly accidental, however, and is due to the nonlinear bias varying above 
and below the linear one depending on the length scale, leading to partial cancellation of the differences when the correlation 
function is integrated over the length scales. For small observational bandwidths, Afobs, there is less cancellation and the 
differences in the two predictions grow at all beam sizes, ASbeam. Similarly, if ASbeam is small, there is little cancellation, 
and the discrepancies are larger at all bandwidths. In these cases the linear bias gives an overestimate of the 21-cm emission 
fluctuations from minihalos at the low end of the redshift range, and an underestimate at high-z, even at large values of Ai^obs- 
However, when both the beam size and the frequency bandwidth are large the differences in the two approaches at small scales 
are diluted, and the resulting 21-cm emission fluctuations become identical. Finally, we predict that the best observational 
frequency bandwidth for improving the chances for detection is Ai/obs ~ 2 MHz. Thus it may be through a such frequency 
window that astronomers get their flrst glimpses of the cosmological Dark Ages. 
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